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We propose global hybrid approximations of the exchange-correlation (XC) energy functional 
which reproduce well the modified fourth-order gradient expansion of the exchange energy in the 
semiclassical limit of many-electron neutral atoms and recover the full local density approximation 
(LDA) linear response. These XC functionals represent the hybrid versions of the APBE functional 
[Plrys. Rev. Lett. 106, 186406, (2011)] yet employing an additional correlation functional which 
uses the localization concept of the correlation energy density to improve the compatibility with 
the Hartree-Fock exchange as well as the coupling-constant-resolved XC potential energy. Broad 
energetical and structural testings, including thermochemistry and geometry, transition metal com¬ 
plexes, non-covalent interactions, gold clusters and small gold-molecule interfaces, as well as an 
analysis of the hybrid parameters, show that our construction is quite robust. In particular, our 
testing shows that the resulting hybrid, including 20% of Hartree-Fock exchange and named hAPBE, 
performs remarkably well for a broad palette of systems and properties, being generally better than 
popular hybrids (PBE0 and B3LYP). Semi-empirical dispersion corrections are also provided. 


I. INTRODUCTION 

Density functional theory (DFT)i^ is nowadays one of 
the most popular methods in electronic structure calcu¬ 
lations. DFT is an exact theory, but its practical im¬ 
plementation requires an approximation for the so-called 
exchange-correlation (XC) energy functional, which de¬ 
scribes the quantum effects of the electron-electron in¬ 
teraction. Many approximated expressions, having dif¬ 
ferent levels of complexity/sophistication, exist for this 
functional^. They can be roughly classified in two broad 
classes: i) local and semilocal functionals, using as in¬ 
put information the electron density, its derivatives, and 
the kinetic energy density; ii) fully non-local function¬ 
als, which additionally use as input non-local quanti¬ 
ties computed from the Kohn-Sham orbitals (e.g. ex¬ 
act exchange). The former ones, which include local 
density, generalized gradient and meta generalized gradi¬ 
ent approximations (LDA, GGAs, and meta-GGAs, re¬ 
spectively), are fast and sufficiently accurate for many 
purposes^—, but fail for a number of important proper¬ 
ties such as density distributions or barriers heights^ - —. 
The second ones, which include global, range-separated, 
local and double hybrid functionals 22 ^ - —, improve con¬ 
siderably the results, but they are much more computer¬ 
time demanding, in particular when local or double hy¬ 
brids are used. 


Among functionals of the second class, global hybrids 
provide the best compromise between accuracy and ef¬ 
ficiency. They are characterized by an exchange energy 
contribution which is a mixing of exact exchange and a 
local or semilocal contribution. In global hybrids, a frac¬ 
tion of Hartree-Fock exchange energy is directly com¬ 
bined with the complementary fraction of an exchange 
energy functional. 

The great majority of the hybrid functionals are 
parameterized: the constant determining the relative 
weight of the exact and semilocal exchange and (in most 
cases) some other parameters are fixed by fitting some 
training datasets. Furthermore, fitted parameters can 
also be contained in the semilocal functional used in the 
hybrid construction. A well-known example of parame¬ 
terized hybrid functional is B3LYP—. 

On the other hand, there are few examples of non- 
parameterized hybrid functionals. In the case of global 
hybrids (the one of interest in the present paper) 
the parameter determining the ratio between exact 
and semilocal exchange can be chosen by theoretical 
argument s 33 : 36 . Then, constructing the hybrid on the 
basis of a non-parameterized semilocal functional, one 
obtains a parameter-free (or non-empirical) hybrid func¬ 
tional. PBEO^ - — and, more recently, PBEO-1/3 36 ’ 37 
have been obtained in such a way. Extensive calcu¬ 
lations of molecular properties have shown that non- 
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parameterized hybrids are not only more satisfactory 
than the parameterized ones from the theoretical point of 
view, but they are also competitive in terms of accuracy, 
in particular when systems or properties that do not be¬ 
long to the training sets of the parameterized hybrids are 
considered. 

In the present paper, we construct and test parameter- 
free hybrids based on the APBE GGA functional^. This 
is a non-parameterized GGA constructed from the semi- 
classical atom model and has shown high accuracy for a 
wide range of molecular properties^. Moreover, in our 
development of the hybrid form, we deviate from the con¬ 
ventional scheme and explore the possibility of combining 
the exact exchange with an additional correlation func¬ 
tional. This procedure is motivated, within the adiabatic 
connection schem o 26 ' 33 ! 36 , by the need to correct the cor¬ 
relation contribution associated to the exact exchange 
fraction, when the coupling constant is approaching one. 
A similar motivation stands at the base of double hy¬ 
brids, where the correlation is however partly treated at 
the MP2 level. 

The paper is organized as follows: in the next section 
we present the theory at the base of our construction 
of the hybrid functionals and discuss it in the context 
of the adiabatic connection. A possible spin-dependent 
correction for the correlation is discussed as well. In the 
following section we show the performance of our hybrid 
functionals in comparison with other relevant ones, using 
a broad set of molecular tests. The results are discussed 
also in terms of the role of the parameters defining the 
hybrids. Finally, we provide a brief study on the use of 
a semiempirical dispersion correction in conjuction with 
the functionals presented in this work. 

II. THEORY 

A. Construction of the hybrid functional 

According to the adiabatic connection the exchange- 
correlation (XC) functional is given by the coupling- 
constant integration formula 

E xc = [ W XCi \dX , (1) 

Jo 

where W xc , a is an appropriate coupling-constant-resolved 
XC potential energy. One popular choice for W xc> \ to 
derive global hybrid functional is the ansatz proposed by 
Perdew, Ernzerhof, and Burke^: 

w xc , A = W?£ A + {E* F - E FFA ) (1 - A)" -1 , (2) 

where DFA stands for some local or semilocal density 
functional approximation and E^ F is the Hartree-Fock 
exchange computed with Kohn-Sham orbitals. This for¬ 
mula correctly reverts to exact-exchange-only when A = 
0 (i.e. for the non-interacting system) and assumes that 
W xc .i is well approximated by DF A 33 ' 64 . We recall^ 


that W FF X A = E ffa + U FFA { A), where U FFA { A) = 

^{A 2 £? FA [Pt,a-i,Pi,a-i]}, with Per, 7 (r) = 7 3 Pcr(7 r ) 
being the uniform scaled spin-densities. Thus, U FFA ( A) 
is a density functional approximation for the correlation 
potential {U? FA { 0) = 0 and U FFA (l) = E FFA -T FFA ), 
and Eq. © describes, at small non-vanishing values of 
the coupling constant A, a mixing of the DFT exchange 
with the exact exchange, but with full DFT correlation 
as given by DFA. Using Eq. © in Eq. © we have: 

^ = lEr + (l-IW M + £? M P) 

n \ n J 

which shows that always the same E FFA is used in com¬ 
bination with different value of n. This situation ap¬ 
pears not to be optimal, since in general the DFA cor¬ 
relation contribution is designed to work well together 
with its DFA exchange counterpart but not with the HF 
exchange. Thus, it may be appropriate to generalize Eq. 
©• Hence, we propose the generalized ansatz 

w*c,a = 

+ (E FFA2 -E ffa1 ) A m -\ (4) 

so that after A integration we obtain 

E„ = IfiT + (l - I) Ef™ + (5) 

+-E? FA 2 +(l--) E FFAl . 
m \ m J 

where the integer parameter m controls how the correla¬ 
tion functionals E FFAi and E FFA2 are mixed. We note 
that Eq. © is exact at A = 0, as Eq. m, while at A = 1 
we have W XCi i = W F C FA1 + E FFA2 — E FFA1 , differently 
from Eq. ©. As shown in the Appendix, Eq. © can 
yield an improved coupling-constant-resolved XC poten¬ 
tial for atoms. In this work, we will choose for DFA1 the 
APBE GGA functional^, which is one of the most ac¬ 
curate non-empirical GGAs with broad applicability^^. 
For DFA2, we select the PBEloc correlation functional^, 
which has been proved to work quite well together with 
the Hartree-Fock exchange, even if other possibilities 
could also be considere d 6 ' 67 ' 68 . Note that Eq. © re¬ 
sembles the general expression of two-parameter double 
hybrid functionals^, if E FFA2 is replaced by the MP2 
correlation energy and the constraint on n and m to be 
integer numbers is relaxed. 

In order to use only one parameter, we need a relation 
between n and to. We note that the APBE functional was 
constructed in such a way that it recovers the accurate 
LDA linear response behavio r 69 ' 70 . Similarly the PBE¬ 
loc functional was designed to respect the same property 
with respect to exact exchange. Thus, it appears a nat¬ 
ural choice to force the full hybrid functional to respect 
the accurate LDA linear response behavior. To this end 
we need to put to = n in Eq. ©. The final formula for 
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the hybrid XC functional is therefore 

E xc = - (E bf + e FBEI oc ) + (1 - -) E^ PBE . (6) 
n v \ n J 

Equation ([6]) shows that the resulting XC energy is just 
a linear mixing of the original APBE XC functional with 
the E x F +E PBEloc functional: thus a good accuracy can 
be expected only if PBEloc works well in combination 
with HF, as it is indeed the case^. 

We recall that the relevance of the LDA linear response 
behavior has been already discussed for the semilo¬ 
cal level of the theor y---^ — and even for hybrid 
functionals^. However, in Ref. [73j only the linear re¬ 
sponse behavior of the semilocal fraction of the functional 
is considered. We will discuss this item even further in 
the final part of this paper (Section IIV A[) . where the 
condition n = m will be released. 

Moreover, we recall that E£ pbe was constructed from 
the semiclassical theory of large neutral atoms, and 
recovers the modified second-order gradient expansion 
(MGE2)£2. It was proved that MGE2 can reproduce well 
the fourth-order gradient expansion (GE4) 62 i 74 i 75 . Thus, 
because E BF already recovers the GE4, the whole ex¬ 
change part of the hybrid given in Eq. © reproduces to a 
good accuracy the GE4 in the slowly-varying atomic den¬ 
sity limit. This feature however, can be met exactly only 
at the meta-GGA hybrid level (e.g. the TPSSh hybrid^). 

The value of the remaining parameter n can be 
strongly circumscribed by theoretical considerations 
based on perturbation theor y 33 : 36 but not univocally 
fixed. The more convenient choice depends on the semilo¬ 
cal functional(s) used in the hybrid construction and is 
ultimately determined by the result accuracy. In this 
work we consider n = 4 or 5. The corresponding func¬ 
tionals are assessed in next sections. These functionals 
have the same exact exchange content as the popular 
PBEO^c— and B3LYP—hybrid functionals, respec¬ 
tively. We note that additional values of n can be con¬ 
sidered. In particular, the value n = 3, which can also 
be obtained from perturbation theory2&, was shown to 
yield good results for organic chemistry^! and excited 
states, in particular when charge transfer 77 and Rydberg 
transitions^ are considered. However, such higher frac¬ 
tions of exact exchange do not provide a global result 
improvement for the datasets considered in the present 
paper, thus they will not be discussed any longer in this 
work. 


B. Spin-dependent correction for the PBEloc 
functional 

The good compatibility of the PBEloc functional with 
the Hartree-Fock exchange is based on a localization 
para digm of the correlation energy density ^ 10 ' 66 . In 
Ref. ^7| it has been shown that this localization, and 
consequently the compatibility with exact exchange, can 
be further enhanced, preserving the exact properties of 





FIG. 1: E c (£) versus £ of the one-electron Gaussian and Hy- 
drogenic densities with uniform spin-polarization £, for PBE¬ 
loc, zvPBEloc, and ideal E c (z) of Eq. ©. 

PBEloc, by a slight modification of the original correla¬ 
tion functional to 

E zvpbeioc = J p [ Y )f a u [Q(r) iV {r))e PBEloc {r)dY , (7) 

where p is the electron density, e PBEloc jg the PBEloc 
correlation energy per particle, and 

/a, w (C,t;) = e-“ u3 ^“ J (8) 

is a spin-dependent correction factor with v = 
\S7n\/{2k v p) being a reduced gradient suitable for den¬ 
sity variations in valence and bonding regions^, k v = 
2(3/(47r 4 )) 1 / 18 p 1 / 9 , and £ being the relative spin polar¬ 
ization. 

To fix the parameters a and ui in Eq. © we used the 
uniformly spin-polarized Gaussian (G) and hydrogenic 
(H) one-electron density models n-j- = ug,h and nj. = 

~Y^n G ,H, with n H {r) = and n G (r) = respec¬ 
tively, requiring an improved spin-behavior for £ > 0.3. 
We recall that the ideal spin-behavior for these model 
systems is^ 

E c (0 = E pbei ° c (( = 0 )g(C). (9) 

with <?(£) = 1 — £ 2 . We finally obtain a = 0.5 and 
u> = 2.0. As shown in Table[j]the so obtained correlation 
zvPBEloc functional displays an improved compatibility 
with Hartree-Fock exchange, with respect to the original 
PBEloc, for tests on small molecules where dynamical 
correlation is relevant (see Refs. I67ll96l for more details). 
Moreover, Fig. [L] shows that zvPBEloc provides a bet¬ 
ter behavior with respect to PBEloc, yet preserving the 
PBEloc shape, for both H and G densities. 

To take advantage of the improved features of the zvP¬ 
BEloc functional we will thus consider also global hybrids 
of the form 

E xc = - (E bf + E z c vPBEloc ) + (l - E x pbe . (10) 
n v \ n J 
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III. COMPUTATIONAL DETAILS 

To assess the performance of the functionals we per¬ 
formed calculations on a large test set (overall more 
than 500 systems and various energetical and structural 
ground-state properties), which was divided into five 
main groups: 

• Thermochemistry. This group considers 

the atomization energies of small molecules 
(AE 6 82 i 83 and W l 55 i 85 i 86 ) and molecules with 
non-single-reference character (W4-MR^), 
proton affinities (PAlS^Sr—) and ionization 
potentials (G21IP— i 86 ' 90 ), different barrier 

heights (BH76 85 i 86 i 91 i 92 ), reaction energies 
(BH76RC—and OMRE^), and both 
barrier heights and reaction energies (K9 83 i 84 ). 

• Organic molecules’ geometry. This group 
comprises tests on bond lengths of hydrogenic 
(MGHBL922) and non-hydrogenic (MGNHBLU^) 
bonds, bond lengths of open-shell organic molecules 
(BLO^l) as well as vibrational frequencies (F3822) 
of small organic molecules. 

• Transition metals. This group of tests in¬ 
cludes atomization energies of small transition 
metal complexes (TMIOA E 95 ! 96 ) and gold clusters 
(AUnA E 96 ! 97 , reaction energies of transition metal 
complexes (TMR E 23 ' 95 ), and bond lengths of tran¬ 
sition metal complexes (TMBL 96 i 98 ) and gold clus¬ 
ters (AuBL6 — > — ). Note that here for the TMIOAE 
test have been used tighter convergence criteria 
than in previous publications^*^. 

• Non-covalent interactions. This group consid¬ 
ers interaction energies of hydrogen-bond (HB6^), 
dipole-dipole (DI622), charge-transfer (CT7—), and 
7r-7r-stacking (pp522) complexes. In addition, in¬ 
teraction energies of dihydrogen-bond complexes 
(DHB23i2£) and of complexes with various char¬ 
acter (S2 2 ' 194 ’-'fi3 .) are taken into account. 

• Other properties. This group considers a miscel¬ 
laneous of other properties and systems including 


isomerization energies of large organic molecules 
(ISOL6i!E), difficult cases for DFT (DC9/12M), 
small gold-organic interfaces (SI1322), dipole mo¬ 
ments of organic molecules (DM25^), and atomic 
energies (AE17—). 


All tests were carried out for the functionals defined 
by Eqs. © and m with n = 4 or 5. For completeness 
also the hybrid functional named APBEO and defined as 
E £ PBE0 = 0.25E bf + 0.75E£ pbe + E^ pbe was consid¬ 
ered (this is in practice obtained from Eq. © changing 
E pbeioc with e apbe and setting n = 4 ). 

For comparison we performed calculations using sev¬ 
eral other GGA and hybrid functionals: the PBE^i and 
APBE22 GGA functionals, as well as the popular hybrid 
XC functionals PBEOH - 35 and B3LYP^rH. Finally, we 
considered also the PBEmol/30 hybrid functional^ which 
was constructed to restore the LDA linear response in the 
semilocal part of the functional by scaling the second- 
order coefficient /3 in the correlation part. Note that, 
because the GGA functional PBEmol 73 is very similar 
to APBE, the PBEmol/30 hybrid functional practically 
differs from Eq. © only for the fact that in the former 
the LDA linear response condition is enforced only in 
the semilocal part, while in the later it is extended to the 
whole functionals (thanks to the inclusion of the PBEloc 
correction). 

All calculations have been performed with the TUR¬ 
BOMOLE program package^ using a def2-TZVPP ba¬ 
sis se t 106,107 . The choice of the basis set was motivated 
by the need to find a best compromise between accuracy 
and computational cost, so that our results can be di¬ 
rectly compared to practical applications. We note also 
that because the same basis set is used for all the func¬ 
tionals a fair assessment of relative performances is pos¬ 
sible. 

To evaluate the performance of the different ap¬ 
proaches we computed, for each group of tests outlined 
above, the average MAE relative to PBEO, which is as¬ 
sumed as a reference. Hence, we considered 


TABLE I: Mean absolute errors (kcal/mol) on AE6 tesl42 
for atomization energies of organic molecules, BH6 tesl4^ for 
barrier heights of organic molecules, and I<9 test^ for kinet¬ 
ics of organic molecules, as resulting from calculations using 
Hartree-Fock exchange and different GGA correlation func¬ 
tionals. 


Test set 

PBE 

LYP 

PBEloc 

zvPBEloc 

AE6 

31.9 

38.2 

24.0 

20.8 

BH6 

5.6 

5.3 

4.4 

4.6 

K9 

5.7 

6.0 

4.7 

4.2 


M 


1 A MAE. 
RMAE = — V 


M MAE: 
1=1 


PBEO 


(ii) 


where the sum runs over all the M tests within a 
group and MAE BBB0 is the MAE of PBEO for the i-th 
test. The RMAE indicates whether any method is bet¬ 
ter (R.MAEcl) or worse (RMAE>1) than PBEO. The 
RMAE provides a fair global assessment of all the re¬ 
sults, however it may tend to overweight tests where both 
MAE* and MAEf BB0 are small and to underweight the 
results of tests where both methods yield quite poor re¬ 
sults. For this reason, in addition to the RMAE, we 
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considered also the weighted MAE 


WMAE 


1 

M 

1 

M 


^ MAE* - MAEf ' BE0 
h <MAEf BB0 ) " ,L 

^ / MAEj \ MAE BBB0 
^ VMAEf BB0 ~~ ) (MAEf BB0 ) ’ 


where (MAE BBB0 ) is the average of the MAEs of all hy¬ 
brid methods for test i. Equation (TH21) shows that in the 
WMAE the quantities composing the RMAE are aver¬ 
aged with weights Wi = MAE BBB0 /(MAE BBB0 ), which 
indicate the relative performance of PBEO (assumed as a 
reference) with respect to what may be expected from hy¬ 
brid approaches. Hence, negative WMAEs denote meth¬ 
ods outperforming PBEO, while positive values denote 
methods performing worse than PBEO. 


IV. RESULTS 

Table M reports the mean absolute error (MAE) for 
each test as obtained from the different functionals. We 
see that the best overall performance is given by the hy¬ 
brid functional defined by Eq. m with n = 5 (named 
hAPBE hereafter), which yields a global RMAE of 0.90 
and a global WMAE of -13%. Notably it is also the 
best non-fitted hybrid functional for three out of five 
groups of tests that we considered (thermochemistry, or¬ 
ganic molecules’ geometry, and transition metals) and it 
is very close to the best results for the remaining two 
groups. Moreover, hAPBE provides results more accu¬ 
rate than the average of the considered hybrids in 23 over 
28 cases. For thermochemistry the hAPBE functional 
yields performance very close to the B3LYP functional 
(RMAE is the same; WMAE is only 2% better in B3LYP) 
which was fitted on some of these properties. However 
B3LYP shows, several limitations for other problems, es¬ 
pecially those related to transition metals and clusters, 
being below the average performance in most tests not 
connected to thermochemistry. These limitations may 
possibly trace back to the fitting parameters as well to 
some limitations of the Lee-Yang-Parr correlation^ in 
the slowly-varying density limit!2£. 

We also note that the hAPBE functional is also supe¬ 
rior to some non-empirical meta-GGAs, such as TPPSHiL 
and BLOC & 10 i 66 which have total WMAEs +16 and +13 
respectively, and simple hybrid meta-GGAs (TPPSh^ 
has WMAE +3). It is still not as accurate as the most 
sophisticated hybrid meta-GGAs. Hence, for example 
the M06 functional^! performs on the test set used in 
this work with a WMAE of -17%. However, we need to 
stress that the M06 functional contains several tens of 
empirical parameters. 

To understand better the results of Tab. El two main 
trends are worth of an analysis. The first one concerns 
the effect of the inclusion of different amounts of Hartree- 
Fock exchange. Information on such a trend can be ob¬ 


tained comparing APBE with APBEO, PBE with PBEO, 
and in general the hybrids with n = 5 and n = 4. 
The analysis of the data in Tab. El shows that in fact 
the inclusion of a fraction of Hartree-Fock exchange can 
benefit the overall performance, but a larger amount of 
Hartree-Fock exchange does not necessarily correspond 
to an improvement of the results. More in detail we 
see that non-covalent complexes are almost insensitive 
to the inclusion of Hartree-Fock exchange (with the par¬ 
tial exception of charge-transfer complexes) whereas or¬ 
ganic molecules and transition-metal complexes have dif¬ 
ferent and quite erratic behaviors. Hence, for organic 
molecules Hartree-Fock exchange provides an improve¬ 
ment for barrier heights and most reaction energies. On 
the other hand, the inclusion of Hartree-Fock exchange 
yields a clear worsening of transition-metal complexes en¬ 
ergies, but a moderate improvement for the bond lengths. 
Thus, a proper balance, taking into account all these ef¬ 
fects, appears hard to find. Nevertheless, it seems that 
a best overall performance may be obtained for a mod¬ 
erate fraction of Hartree-Fock exchange. This result is 
in agreement with the finding of Ref. [23|. On the other 
hand, even within functionals having the same amount of 
Hartree-Fock exchange, significant differences are found. 
In particular, we remark that from the present assess¬ 
ment the hybrids proposed in this work appear to per¬ 
form better than the popular hybrids B3LYP and PBEO. 

The second trend to observe concerns the different pos¬ 
sible choices for the correlation part of the functional. 
In particular, it is interesting to inspect the APBEO, 
PBEmol/30, and the functional defined by Eq. © with 
n = 4 (the spin-dependent correction for PBEloc is dis¬ 
cussed below), to understand the role of the LDA linear 
response constraint and of the PBEloc correlation. In 
agreement with Ref. [73j we find that the imposition of 
the LDA linear response on the semilocal part of the 
functional gives in general an improvement of the perfor¬ 
mance, especially for atomization and reaction energies 
of organic molecules. However, in the case of PBEmol/30 
this seems to occur mostly thanks to an error cancel¬ 
lation effect, as shown by the fact that atomic energies 
are much worse for PBEmol/30 than for APBEO. On the 
other hand, the functional proposed in this work shows a 
similar (but systematically better) performance as that of 
PBEmol/30, with apparently a smaller error cancellation. 

We remark once more that neither APBEO nor 
PBEmol/30 recover the full LDA linear response, as in¬ 
stead do the functionals of the present work. In fact, the 
exchange-correlation second-order gradient coefficients 
are: /i BBB0 = -0.024, ^ BBB0 = -0.034, /r BBBmo ^ 0 = 
0.031, while only for the functionals of Eqs. © and (9), 
H xc = 0. Because in Table II the latter functionals give 
the best global RMAE and global WRMAE, the imposi¬ 
tion of the full LDA linear response for the full global hy¬ 
brid seems to be important. This correlation issue may 
appear, at first sight to have only a minor impact on 
most calculations, because energy differences are often 
considered. However, it may be relevant for those cases 
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TABLE II: Mean absolute errors (MAEs) for different tests and functionals. For each group of tests the average MAE relative 
to PBEO and the weighted MAE (RMAE and WMAE; see text for definitions) are reported in the last lines. For each line the 
hybrid functional performing best is highlighted in bold face; a star is used to denote the best among the functionals based on 
Eqs. © and ©. 


Test 

GGA 

Hybrid n = 

= 5 


Hybrid n = 4 




PBE 

APBE 

B3LYP Eq. © Eq. (fTUll 

PBEO 

APBE0 PBEmol/?0 Eq. © Eq. (USD 

atomization energy (AE6) 

13.4 

Thermochemistry (ki 
7.8 5.5 5.8 

:al/mol) 

* 4.2 

6.3 

9.5 

7.3 

6.8 

4.8 

atomization energy (W4) 

10.7 

8.5 

5.8 

6.4 

*6.0 

6.0 

8.7 

7.6 

7.3 

6.4 

static correlation (W4-MR) 

21.8 

17.4 

8.4 

* 5.9 

* 5.9 

7.1 

10.5 

11.2 

9.1 

7.2 

proton affinities (PA12) 

2.2 

2.8 

2.3 

*3.1 

*3.1 

2.8 

3.4 

3.5 

3.3 

3.3 

ionization pot. (G21IP) 

3.9 

4.0 

3.8 

*4.0 

*4.0 

4.1 

4.2 

4.6 

4.2 

4.2 

barrier heights (BH76) 

9.8 

9.0 

5.0 

5.0 

5.0 

4.5 

3.9 

4.0 

*4.2 

*4.2 

reaction energies (BH76RC) 

4.4 

3.8 

2.6 

2.8 

2.9 

2.7 

2.5 

2.4 

*2.6 

2.7 

reaction energies (OMRE) 

6.7 

7.1 

4.6 

6.5 

*6.4 

9.1 

7.1 

6.0 

8.1 

7.9 

kinetics (K9) 

7.5 

6.6 

3.8 

4.2 

4.2 

3.9 

3.3 

3.1 

*3.7 

*3.7 

RMAE 

1.71 

1.46 

0.96 

1.00 

* 0.96 

1.00 

1.15 

1.09 

1.07 

0.99 

WMAE(%) 

+66 

+42 

-6 

-1 

*-4 

0 

+11 

+6 

+6 

-2 

H bond lengths (MGHBL9) 

Organic molecules’ geometry 
11.4 10.2 3.0 *1.4 

(mA, and 
*1.4 

cm" 1 ) 

2.3 

1.5 

1.0 

1.5 

1.5 

non-H bond lengths (MGNHBL11) 

7.6 

9.2 

7.2 

*7.3 

*7.3 

9.8 

8.8 

8.7 

10.7 

10.6 

open-shell molecules (BL9) 

15.0 

14.7 

11.9 

* 11.6 

* 11.6 

12.7 

12.2 

12.3 

12.5 

12.5 

vibrations (F38) 

58.4 

55.0 

37.1 

*44.2 

*44.2 

53.3 

52.0 

52.1 

58.3 

58.3 

RMAE 

2.00 

1.88 

0.91 

* 0.77 

* 0.77 

1.00 

0.86 

0.82 

0.96 

0.96 

WMAE(%) 

+135 

+119 

-7 

*-28 

*-28 

0 

-17 

-23 

-7 

-7 

Transition metal complexes 
atomization energy (TM10AE) 13.0 11.1 

(kcal/mol 

13.4 

[for AUnAE kcal/(mol-atoms)] and 
12.1 * 11.0 14.4 15.8 

mA) 

14.0 

14.1 

11.4 

reaction energies (TMRE) 

3.7 

3.1 

10.6 

9.6 

* 8.4 

9.0 

9.5 

9.5 

11.5 

9.2 

gold clusters atomiz. (AUnAE) 

0.60 

1.8 

5.9 

* 3.8 

* 3.8 

4.0 

5.1 

5.0 

4.3 

4.3 

bond lengths (TMBL) 

13.5 

13.0 

18.5 

* 18.3 

* 18.3 

21.2 

20.6 

20.5 

21.8 

21.8 

gold clusters bonds (AuBL6) 

56.5 

58.7 

94.2 

36.6 

36.5 

41.9 

56.0 

58.9 

32.5 

* 32.3 

RMAE 

0.69 

0.72 

1.34 

0.92 

* 0.88 

1.00 

1.15 

1.13 

1.03 

0.94 

WMAE(%) 

-30 

-28 

+29 

-8 

*-12 

0 

+13 

+11 

+3 

-6 

hydrogen bonds (HB6) 

0.38 

Non-covalent interactions (kcal/mc 
0.32 0.55 * 0.36 * 0.36 

'1) 

0.52 

0.36 

0.45 

0.38 

0.38 

dipole-dipole (DI6) 

0.38 

0.39 

0.88 

*0.38 

*0.38 

0.36 

0.38 

0.32 

*0.38 

*0.38 

dihydrogen bonds (DHB23) 

0.98 

0.76 

0.34 

*0.66 

*0.66 

0.75 

0.56 

0.79 

*0.66 

*0.66 

charge-transfer (CT7) 

2.7 

2.4 

0.79 

1.3 

1.3 

1.3 

1.0 

1.2 

*1.2 

*1.2 

7r-7r stacking (pp5) 

2.1 

2.2 

3.2 

*2.2 

*2.2 

2.1 

2.2 

2.0 

*2.2 

*2.2 

various non-covalent (S22) 

2.3 

2.7 

3.5 

*2.4 

*2.4 

2.2 

2.4 

2.2 

*2.4 

*2.4 

RMAE 

1.21 

1.13 

1.27 

0.97 

0.97 

1.00 

0.91 

0.95 

*0.95 

*0.95 

WMAE(%) 

+23 

+15 

+22 

-5 

-5 

0 

*-12 

-6 

*-7 

*-7 

isomerization (ISOL6) 

2.2 

Other (kcal/mol; 
2.4 2.7 

Debye/10 for DM25) 

* 1.4 * 1.4 1.5 

1.4 

1.4 

* 1.4 

* 1.4 

difficult cases (DC9.12) 

40.8 

29.4 

25.0 

19.7 

19.6 

17.7 

22.9 

20.2 

18.2 

*17.9 

small interfaces (SI12) 

3.7 

5.9 

11.0 

* 7.0 

7.1 

7.2 

8.9 

9.3 

7.8 

7.9 

dipole moments (DM25) 

3.6 

3.6 

2.9 

2.7 

2.7 

2.7 

2.7 

2.7 

* 2.6 

* 2.6 

atomic energies (AE17) 

51.6 

22.2 

7.6 

16.9 

16.6 

42.2 

10.2 

54.5 

15.7 

*15.3 

RMAE 

1.36 

1.19 

1.20 

0.88 

0.88 

1.00 

0.94 

1.13 

0.88 

* 0.87 

WMAE(%) 

+38 

+9 

+2 

-22 

-23 

0 

-21 

+17 

*-24 

*-24 

Global RMAE 

1.41 

1.27 

1.12 

0.92 

* 0.90 

1.00 

1.02 

1.03 

0.99 

0.94 

Global WMAE(%) 

+44 

+28 

+5 

-11 

*-13 

0 

-4 

+2 

-5 

-9 
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TMRE 




i 

■ 


0.88997 

0.88630 

0.88263 

0.87896 

0.87530 

0.87163 

0.86796 

0.86429 


0.85329 

0.84962 

0.84595 

0.84228 

0.83861 


FIG. 2: Mean absolute errors (kcal/mol) for several energy tests and different values of the parameters in the hybrid functional 
of Eq. (1131) . The scaled average (see text for the definition) of the various tests is also reported (right bottom panel). In each 
panel the white boxes denote the positions of the functionals defined by Eq. m with n = 4 and n = 5 as well as APBEO. 


where heterogeneous systems are involved, such as for 
the atomization of metal clusters or the description of 
metal-organic interfaces, because in these cases the cor¬ 
relation effects of the different parts are likely to be very 
different and will not cancel properly. This observation 
is supported by the trends registered in Tab. [TT] for the 
AUnAE and SI12 tests, although the small dimensions 
of the systems considered in those tests prevent a good 
analysis of this effect. 

Finally, we consider a comparison between hybrid func¬ 
tionals including the simple PBEloc correlation and those 
implementing the spin-dependent correction via zvPBE- 
loc. The data of Table mi show the latter functionals 


display a non-negligible improvement with respect to the 
former ones: the improvement in WMAE is 4% (2%) for 
n = 4 (n = 5). This fact can be traced back to the 
higher compatibility with Hartree-Fock exchange of zvP- 
BEloc with respect to PBEloc. In fact, the more signif¬ 
icant improvement is registered for n = 4, while smaller 
advantages are observed for n = 5. 

Also interesting to note is the fact that the inclusion of 
the spin-dependent correction brings in all cases either an 
improvement of the results or leaves them practically un¬ 
changed (a very small worsening is observed only in few 
cases). We recall, in addition, that for spin-compensated 
cases (e.g. all the non-covalent tests) the spin-correction 
































has no effect, by construction. 


A. Analysis of hybrid parameters 

To gain more insight into the performance of the hybrid 
functionals we consider in this section the general expres¬ 
sion in Eq. (fojl, with DFA1=APBE and DFA2=PBEloc, 
i.e. 

E xc = aE BF +{l-a)E^ PBE +PE PBEloc +{l-p)E^ PBE , 

(13) 

and we preform a scan of the values of the two parameters 
a and /3 to investigate the dependence of the results on 
the fraction of Hartree-Fock exchange and the correlation 
contribution. To this end we consider a minimal set of 
tests composed of the energy tests AEG (atomization en¬ 
ergies of small molecules), K9 (barrier heights and reac¬ 
tion energies of organic molecules), and TMRE (reaction 
energies of transition-metal complexes). These tests are 
in fact representative of the most important systems and 
properties as well as of the different trends observable in 
Table [XT] To evaluate a global performance of the energy 
tests, which display rather different MAEs, we consider 
an average of the three tests after rescaling the result of 
each test as MAE—>MAE/(1+MAE). Note that the scal¬ 
ing has the two-fold purpose of making all the quantities 
comparable, so that a simple average is meaningful, and 
to enhance the differences between small values (which 
are the most interesting here). 

The results of the scan for the energy tests are reported 
in Fig. [2] The plots show that indeed, as noticed when 
discussing Tab. El the different tests display different be¬ 
haviors with respect to a and /?. In particular, the bar¬ 
rier heights within the K9 test require a large fraction of 
Hartree-Fock exchange for an accurate description, while 
the opposite occurs in transition metals. On the other 
hand, atomization energies of organic molecules show an 
even more complex trend, requiring a delicate balance 
between exact exchange and correlation. Thus, overall 
the average performance has a complicated dependence 
on the two parameters. Nevertheless, it is possible to 
identify a clear minimum approximately corresponding 
to a = 0.2 and /? = 0.3, which is close to the definition of 
the functional with n = 5. Moreover, reasonably small 
average errors can be obtained for the family of parame¬ 
ters respecting the relation /3 ps a + 0.1, which is close to 
the condition used in Eq. ©• Such a relation also shows 
that, as discussed in Section [II A1 an accurate correlation 
functional to be used in hybrid functional should depend 
on the HF exchange contribution. 

These results indicate the robustness of our construc¬ 
tion based on the ansatz of Eq. © and the three imposed 
constraints. In particular, it is highlighted the impor¬ 
tance of the satisfaction of the LDA linear response for 
the full XC functional, which is the constraint determin¬ 
ing the condition m = n in the functional of Eq. ©• 

A somehow different situation is found concerning the 
geometry tests (MGHBL9 and MGNHBL11), which are 


0.30 


8 

0.20 



1 16.19 
15.30 
14.41 
13.52 
12.63 
— 11.74 
10.85 
9.96 
9.07 
8.19 
7.30 
6.41 
5.52 

■ 4.63 
3.74 


FIG. 3: Mean absolute errors (mA) for geometry tests as 
obtained for different values of the parameters in the hybrid 
functional of Eq. m • The white boxes denote the positions 
of the functionals defined by Eq. m with n = 4 and n = 5 
as well as APBE0. 


reported in Fig. [3] In this case in fact a large and shallow 
minimum is observed for moderate values of the param¬ 
eter a and (3 values ranging from 0.2 to 1. This indicates 
that reasonably small errors can be achieved by hybrid 
functionals including a sufficiently moderate fraction of 
Hartree-Fock exchange. In fact, for n = 5, the result of 
4.4 mA is reasonably close to the global minimum of 3.7 

mA. 

Note also that the lines /? = 1 in both Figs. [2] and 
[3] which correspond to the hybrids which use the full 
PBEloc correlation, show a modest performance for all 
the tests. Such hybrids not only violates the LDA linear 
response, but also behave in the tail of the density as 
a pure exchange functional (because of the PBEloc con¬ 
struction). We recall that the correlation energy density 
plays an important role for asymptotic propertie s 109 ’ 110 . 
Nevertheless, remarkably both Figs. [5] and [3] show that 
the best results are found for /3 ss 0.2, which reveals 
an important significance of the PBEloc correlation as 
a tool for the construction of hybrid functionals. Thus, 
further development of semilocal correlation functionals 
more compatible with HF exchange, may further improve 
the quality of global hybrids. 


B. Semiempirical dispersion correction 

To complete the construction of the hybrid functionals 
we considered to complement them with a semiempirical 
dispersion correction 111 . In fact, the dispersion inter¬ 
action cannot be described neither by the Hartree-Fock 
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TABLE III: Numerical values of the parameters used in the 
semiempirical D3 dispersion correction of the different func¬ 
tionals. 


Functional 

Sr,6 

ss 

APBE 

1.242 

0.930 

APBEO 

1.259 

0.921 

this work n = 4 

1.266 

0.915 

this work n = 5 

1.258 

0.923 


exchange nor by any semilocal correlation. We imple¬ 
mented the dispersion correction through the DFT-D3 
semiempirical modeli**, fixing the two free parameters 
of the model by fitting to the MAE of the S22 test. The 
resulting parameters for different functionals are listed in 
Table m Note that because the S22 test only includes 
systems without any spin polarization the parameters of 
the D3 correction are the same for functionals defined by 
Eqs. (0 and (fTOl) . 

The D3 dispersion correction is found to be well com¬ 
patible with all the hybrid functionals, bringing a good 
improvement for dispersion dominated cases (i.e. S22 
and pp5, where van der Waals interactions dominate or 
rather large systems are considered; see Tab. ED. Inter¬ 
estingly the semiempirical dispersion correction also re¬ 
duces the MAE for ISOL6 and metal-organic interfaces 
(SI12). On the other hand, for the first four test in Tab. 
II VI where dispersion is not dominating, the D3 correction 
tends to slightly worsen the results for APBE and its hy¬ 
brid variants. For all other tests (not reported Tab. ITVl) 
where dispersion does not play any relevant role, results 
are correctly left almost unchanged (data non reported). 

V. CONCLUSIONS 

We have developed the global hybrids of the APBE 
exchange-correlation GGA functional^. The construc¬ 
tion of these hybrids is based on the use of the PBE- 
loc or zvPBEloc correlation functionals, to accompany 
the fraction of exact exchange, and on the recovery of 
the accurate LDA linear response. By a broad ener¬ 
getical and structural testing, including thermochem¬ 
istry, organic geometry, atomization energies, reaction 
energies, and bond lengths of transition metal com¬ 
plexes, non-covalent interactions, isomerization, gold- 
molecules hybrid-interfaces, and dipole moments of or¬ 
ganic molecules, we have shown that the best total per¬ 
formance is obtained by considering the hybrid using the 
zvPBEloc correlation (i.e Eq. (1101) ) and n = 5. We name 
this functional hAPBE. The hAPBE hybrid functional, 
shows an almost systematic improvement over the popu¬ 
lar PBEO and B3LYP hybrids, as well as over the recently 
proposed PBEmol/30 functional. In fact, the hAPBE 
hybrid functional has a good accuracy for all the tests, 
showing a broad applicability, in contrast to the B3LYP 


functional which is modest for the transition metal com¬ 
plexes. Use of a semiempirical dispersion correction can 
bring further accuracy for problems where the dispersion 
interaction is especially relevant. Furthermore, the anal¬ 
ysis of the hybrid parameters, summarized in Figs. [2] 
and [31 has shown that the LDA linear response is a pow¬ 
erful constraint, and that the PBEloc correlation plays 
a significant role in the functional performance. Thus 
the hAPBE functional can be regarded as an interest¬ 
ing tool for quantum chemistry applications. Whereas, 
the construction presented in this work, using a different 
correlation contribution to accompany the Hartree-Fock 
exchange and fulfilling the LDA linear response, appears 
to be an important strategy to develop more accurate 
hybrid functionals in the future. 

In view of further improvements different semilocal 
correlation functionals with good compatibility with the 
Hartree-Fock exchange can be considered (e.g. GAPloo^ 
or revTCA£). Alternatively, the replacement of GGAs 
with meta-GGAs can be considered in Eq. 0- We recall 
in fact that important recent work at the meta-GGA level 
has appeared. Hence, several accurate non-empirical 
(semi-empirical) functionals have been proposed, such 
as revTPSSiH, regTPSSi^, MGGA-MSil, VT{8, 4}i±£, 
and BLO C 8 i 10 i 66 . All of them represent good candidates 
for the construction of an accurate meta-GGA hybrid, 
using the method proposed in the present paper. In 
particular the BLOC functional uses a correlation term 
(TPSSloo^) which has been derived from the PBEloc 
GGA. Thus, this functional appears as the most natural 
choice for future studies. 


VI. APPENDIX 

In this Appendix we analyze the coupling-constant- 
resolved XC potential energy formula in Eqs. m and 
0 for atomic systems. As a reference we consider the 
coupling-constant-resolved XC potential energy of the 
interaction-strength interpolation (ISI) model a 

high-level method constructed from exact conditions^. 

The differences between the various approximate W XCt \ 
and the reference, as functions of A, for several neutral 
atoms, are reported in [U In addition we consider also 
the Ne 8+ ion as an example for the high-density limit 
case. For all atoms we report also the exact behavior 
Wax ,a —t E x + 2A E^ L2 , where Ef L2 is the second-order 
Gorling-Levy perturbation theory correlation (GL2) 80 i 81 , 
for small values of A. Note that for Ne 8+ this behav¬ 
ior becomes almost exact over the whole range of A val¬ 
ues. The data in [4] show that the PBEO ansatz is cor¬ 
rect only at A = 0 and in the close proximity of this 
point (by construction), whereas it shows a significant 
inaccuracy for larger A values. The results are much im¬ 
proved for all atoms considered when the semilocal DFT 
approximation is changed from PBE to APBE (APBEO 
curve). In this case in fact a closer behavior to the ref¬ 
erence is obtained over the whole range of A values. By 
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TABLE IV: Mean absolute errors (MAEs) on different tests for the dispersion corrected (D3) hybrid functionals. 



GGA 

Hybrid 

n = 5 


Hybrid n = 4 


Test 

APBE-D3 

Eq. 0-D3 

Eq. (fTUl)-D3 

APBE0-D3 

Eq. @-D3 

Eq. JTDJ-D3 

hydrogen bonds (HB6) 

0.77 

1.0 

1.0 

1.0 

1.1 

1.1 

dipole-dipole (DIG) 

0.91 

0.73 

0.73 

0.69 

0.69 

0.69 

dihydrogen bonds (DHB23) 

1.3 

1.2 

1.2 

1.0 

1.1 

1.1 

charge-transfer (CT7) 

2.9 

1.8 

1.8 

1.5 

1.6 

1.6 

TT-n stacking (pp5) 

0.2 

0.18 

0.18 

0.20 

0.21 

0.21 

various non-covalent (S22) 

0.5 

0.59 

0.59 

0.57 

0.62 

0.62 

isomerization (ISOL6) 

2.2 

1.0 

1.0 

1.0 

1.0 

1.0 

small interfaces (SI12) 

4.5 

5.8 

5.7 

7.6 

6.6 

6.5 



0e+00 

-5e-02 

o 

-le-01 | ^ 

55 i 

£ 

0e+00 hh U 

22 $ 

£ 


-5e-02 


-le-01 






FIG. 4: Difference between several approximate coupling-constant-resolved XC potential energy formulas and the ISI model for 
the He, Be, Ne atoms and the Ne 8+ ion. For small values of A also the GL2 curve is shown as a reference for the initial slope. 


construction, both W PBE0 and W APBE0 are exact at 
A = 0, while at A = 1, they recover the GGA behav¬ 
iors, W X PBE0 -> W PBE and W APBE0 -4 W APBE &. 
As shown in Fig. [U W APBE significantly improves over 
W PBE at A —> 1, and consequently the APBEO hybrid is 
more realistic. 

Then we consider the ansatz introduced in Eq. 0. 


The curves labelled “ Eq.(A)n = 4” and “Eq.( 4)n = 5” in 
U clearly show that the correlation correction has a non 
negligible effect only for A > 0.6 and that, for all system 
considered, a significant improvement is obtained for A 
close to 1. In particular an almost vanishing error for 
A = 1 is obtained for He and Ne. 


Hohenberg, P.; Kohn, W. Phys. Rev. 1964 , 136 , B864- 


Electronic address: eduardo.fabiano@nano.cnr.it 


i 
























11 


B871. 

2 Kohn, W.; Sham, L. J. Phys. Rev. 1965 , 140, A1133- 
A1138. 

3 Perdew J. P.; Schmidt, K. AIP Conf. Proc. 2001 , 577, 
1 - 20 . 

4 Scuseria, G. E.; Staroverov, V. N. In Theory and Applica¬ 
tions of Computational Chemistry: The First 40 Years (A 
Volume of Technical and Historical Perspectives); Dyk- 
stra, C. E.; Frenking, G.; Kim, K. S.; Scuseria, G. E., 
Eds.; Elsevier: Amsterdam, 2005; Chapter 24, pp 669- 
724. 

0 Tognetti, V.; Cortona, P.; Adamo, C. J. Chem. Phys. 
2008 , 128, 034101. 

6 Tognetti, V.; Cortona, P.; Adamo, C. Chem. Phys. Lett. 
2008 , 460, 536-539. 

' Bremond, E.; Pilard, D.; Ciofini, I.; Chermette, H.; 
Adamo, C.; Cortona, P. Theor. Chem. Ace. 2012, 131, 
1184. 

8 Constantin, L. A.; Fabiano, E.; Della Sala, F. J. Chem. 
Theory Comput. 2013 , 9, 2256-2263. 

9 Constantin, L. A.; Chiodo, L.; Fabiano, E.; Bodrenko, I.; 
Della Sala, F. Phys. Rev. B 2011 , 84, 045126. 

10 Constantin, L. A.; Fabiano, E.; Della Sala, F. Phys. Rev. 
B 2013 , 88, 125112. 

11 Armiento R.; Mattsson, A. E. Phys. Rev. B 2005, 72, 
085108. 

12 Mattsson, A. E.; Armiento, R. Phys. Rev. B 2009 , 79, 
155101. 

13 Pittalis, S.; Rasanen, E.; Vilhena, J. G.; Marques, M. A. 
L. Phys. Rev. A 2009 , 79, 012503. 

14 del Campo, J. M.; Gazquez, J. L.; Trickey, S.; Vela, A. 
Chem. Phys. Lett. 2012 , 543, 179-183. 

15 Krieger, J. B.; Chen, .].; Kurth, S. In Density Functional 
Theory and its Applications to Materials; van Doren, V.; 
van Alsenoy, C. P. G., Eds.; Plenum Press: New York, 
2000 . 

16 Zhao Y.; Truhlar, D. G. J. Chem. Phys. 2006 , 125, 
194101. 

17 Sun, J.; Xiao, B.; Ruzsinszky, A. J. Chem. Phys. 2012 , 
137, 051101. 

18 Chiodo, L.; Constantin, L. A.; Fabiano, E.; Della Sala, F. 
Phys. Rev. Lett. 2012, 108, 126402. 

19 Tognetti, V.; Cortona, P.; Adamo, C. Theor. Chem. Ac¬ 
count 2009 , 122, 257. 

20 Tognetti, V.; Adamo, C.; Cortona, P. Interdiscip. Sci. 
Comput. Life Sci. 2010, 2, 163-168. 

21 Janesko, B. G.; Scuseria, G. E. J. Chem. Phys. 2008, 128, 
244112. 

22 Cohen, A. .].; Mori-Sanchez, P.; Yang, W. Chem. Rev. 

2012 , 112, 289-320. 

23 Fabiano, E.; Constantin, L. A.; Della Sala, F. Int. J. 
Quantum Chem. 2013 , 113, 673-682. Fabiano, E.; Con¬ 
stantin, L. A.; Della Sala, F. Int. J. Quantum Chem. 

2013 , 113, 1600. 

24 Grabowski, I.; Teale, A. M.; Smiga, S.; Bartlett, R. J. J. 
Chem. Phys. 2011, 135, 114111. 

25 Grabowski, I.; Teale, A. M.; Fabiano, E.; Smiga, S.; Buk- 
sztel, A.; Della Sala, F. Mol. Phys. 2014 , 112, 700-710. 

26 Becke, A. D. J. Chem. Phys. 1993 , 98, 1372. 

2 ‘ Becke, A. D. J. Chem. Phys. 1996 , 104, 1040. 

28 Savin, A. In Recent Developments and Applications of 
Modern Density Functional Theory, Seminario, J. M. Ed.; 
Elsevier: Amsterdam, 1996; pp. 327-357. 

29 Iikura, H.; Tsuneda, T.; Yanai, T.; Hirao, K. J. Chem. 


Phys. 2001 , 115, 3540. 

30 Vydrov, O. A.; Scuseria, G. E. J. Chem. Phys. 2006 , 125, 
234109. 

31 Zhao Y.; Truhlar, D. G. Theor. Chem. Acc. 2008 , 120, 
215-241. 

32 Zhao, Y.; Schultz, N. E.; Truhlar, D. G. J. Chem. Phys. 
2005 , 123, 161103. 

33 Perdew, J. P.; Ernzerhof, M.; Burke, K. J. Chem. Phys. 
1996 , 105, 9982. 

34 Adamo, C.; Barone, V. J. Chem. Phys. 1999 , 110, 6158. 

35 Ernzerhof M.; Scuseria, G. E. J. Chem. Phys. 1999 , 110, 
5029. 

36 Cortona, P. J. Chem. Phys. 2012 , 136, 086101. 

37 Guido, C. A.; Bremond, E.; Adamo, C.; Cortona, P. J. 
Chem. Phys. 2013 , 138, 021104. 

38 Hermet, J.; Cortona, P.; Adamo, C. Chem. Phys. Lett. 
2012 , 519-520, 145-149. 

39 Haunschild, R.; Janesko, B. G.; Scuseria, G. E. J. Chem. 
Phys. 2009 , 131, 154112. 

40 Haunschild, R.; Perdew, J. P.; Scuseria, G. E. J. Chem. 
Phys. 2012, 137, 224104. 

41 Arbuznikov, A. V.; Kaupp, M. Chem. Phys. Lett. 2007 , 
440 , 160-168. 

42 Arbuznikov, A. V.; Kaupp, M. J. Chem. Phys. 2012 , 136, 
014111. 

43 Arbuznikov, A. V.; Kaupp, M. Int. J. Quantum Chem. 
2011 , 111, 2625-2638. 

44 Haunschild, R.; Henderson, T. M.; Jimenez-Hoyos, C. A.; 
Scuseria, G. E. J. Chem. Phys. 2010, 133, 134116. 

45 Haunschild R.; Scuseria, G. E. J. Chem. Phys. 2010, 132, 
224106. 

46 Perdew, J. P.; Staroverov, V. N.; Tao, .].; Scuseria, G. E. 
Phys. Rev. A 20 08 , 78, 052513. 

47 Odashima, M. M.; Capelle, K. Phys. Rev. A 2009 , 79, 
062515. 

48 Haunschild, R.; Odashima, M. M.; Scuseria, G. E.; 
Perdew, J. P.; Capelle, K. J. Chem. Phys. 2012 , 136, 
184102. 

49 Huang, Y.-W.; Lee, S.-L. Chem. Phys. Lett. 2010 , 492, 
98-102. 

50 Dickson, R. M.; Becke, A. D. J. Chem. Phys. 2005 , 123, 
111101 . 

51 Johnson, E. R.; Becke, A. D. J. Chem. Phys. 2005 , 123, 
024101. 

52 Becke, A. D. J. Chem. Phys. 2005 , 122, 064101. 

53 Johnson, E. R.; Becke, A. D.; Sherrill, C. D.; DiLabio, G. 
A. J. Chem. Phys. 2009 , 131, 034111. 

54 Mori-Sanchez, P.; Cohen, A. .].; Yang W. J. Chem. Phys. 
20 06 , 124, 091102. 

55 Karton, A.; Tarnopolsky, A.; Lamere, J. F.; Schatz, G. C.; 
Martin, J. M. J. Phys. Chem. A 2008 , 112, 12868-12886. 

56 Schwabe, T .; Grimme, S. Phys. Chem. Chem. Phys. 2007 , 
9, 3397-3406. 

57 Grimme S.; Neese, F. J. Chem. Phys. 2007 , 127, 154116. 

58 Becke, A. D. Phys. Rev. A 1998 , 38, 3098. 

59 Lee, C.; Yang, W.; Parr, R. G. Phys. Rev. B 1988 , 37, 
785. 

60 Becke, A. D. J. Chem. Phys. 1993 , 98, 5648. 

61 Stephens, P. .].; Devlin, J. F.; Chabalowski, C. F.; Frisch, 
M. J. J. Phys. Chem. 1994 , 98, 11623. 

62 Constantin, L. A.; Fabiano, E.; Laricchia, S.; Della Sala, 
F. Phys. Rev. Lett. 2011, 106, 186406. 

63 Fabiano, E.; Constantin, L. A.; Della Sala, F. J. Chem. 
Theory Comput. 2011 , 7, 3548-3559. 



12 


64 Arbuznikov, A. V.; Kaupp, M. J. Chem. Phys. 2008, 128, 
214107. 

65 Sharkas, K.; Toulose, J.; Savin, A. J. Chem. Phys. 2011, 
134, 064113. 

66 Constantin, L. A.; Fabiano, E.; Della Sala, F. Phys. Rev. 
B 2012, 86, 035130. 

67 Constantin, L. A.; Fabiano, E.; Della Sala, F. J. Chem. 
Phys. 2012 137, 194105. 

68 Fabiano, E.; Trevisanutto, P. E.; Terentjevs, A.; Con¬ 
stantin, L. A. J. Chem. Theory Comput. 2014, 10, 2016- 
2026. 

69 Moroni, S.; Ceperley, D. M.; Senatore, G. Phys. Rev. Lett. 
1995 , 75, 689. 

70 Ortiz, G. Phys. Rev. B 1992, 45, 11328. 

71 Perdew, J. P.; Burke, K.; Ernzerhof, M. Phys. Rev. Lett. 
1996, 77, 3865. 

72 Cancio, A. C.; Wagner, C. E.; Wood, S. A. Int. J. Quan¬ 
tum Chem. 2012, 112, 3796-3806. 

' 3 del Campo, J. M.; Gazquez, J. L.; Trickey, S. B.; Vela, A. 
J. Chem. Phys. 2012, 136, 104108. 

74 Elliott, P.; Burke, K. Can. J. Chem. 2009, 87, 1485-1491. 

75 Lee, D.; Constantin, L. A.; Perdew, J. P.; Burke, K. J. 
Chem. Phys. 2009, 130, 034107. 

76 Staroverov, V. E.; Scuseria, G. E.; Tao, J.; Perdew, J. P. 
J. Chem. Phys. 2003, 119, 12129. 

77 Guido, C. A.; Cortona, P.; Mennucci, B.; Adamo, C. J. 
Chem. Theory Comput. 2013, 9, 3118-3126. 

78 Guido, C. A.; Cortona, P.; Adamo, C. J. Chem. Phys. 
2014, 140, 104101. 

79 Seidl, M.; Perdew, J. P.; Kurth, S. Phys. Rev. Lett. 2000, 
84, 5070. 

80 Gorling, A.; Levy, M. Phys Rev B 1993, /7, 13105. 

81 Seidl, M.; Perdew, J. P.; Levy, M. Phys. Rev. A 1999, 59, 
51. 

82 Lynch, B. J.; Truhlar, D. G. J. Phys. Chem. A 2003, 107, 
8996-8999. 

83 Haunschild, R.; Klopper, W. Theor. Chem. Ace. 2012, 
131, 1112. 

84 Lynch, B. J.; Truhlar, D. G. J. Phys. Chem. A 2003, 107, 
r 3898-3906. 

85 Goerigk, L.; Grimme, S. J. Chem. Theory Comput. 2010, 
6, 107-126. 

86 Goerigk, L.; Grimme, S. Phys. Chem. Chem. Phys. 2011, 
13, 6670-6688. 

8 ' Parthiban, S.; Martin, J. M. L. J. Chem. Phys. 2001, 
114, 6014. 

88 Martin, .1. M. L.; de Oliveira, G. J. Chem. Phys. 1999, 
111, 1843. 

89 Zhao, Y.; Truhlar, D. G. J. Phys. Chem. A 2006, 110, 
10478-10486. 

90 Curtiss, L. A.; Raghavachari, K.; Trucks, G. W.; Pople, 
J. A. J. Chem. Phys. 1991, 94, 7221. 

91 Zhao, Y.; Lynch, B. J.; Truhlar, D. G. J. Phys. Chem. A 


2004, 108, 2715-2719. 

Zhao, Y.; Gonzalez-Garci'a, N.; Truhlar, D. G. J. Phys. 
Chem. A 2005, 109, 2012-2018. 

Zhao, Y.; Truhlar, D. G. J. Chem. Phys. 2006, 125, 
194101. 

Biczysko, M.; Panek, P.; Scalmani, G.; Bloino, J.; Barone, 
V. J. Chem. Theory Comput. 2010, 6, 2115-2125. 

Furche, F.; Perdew, J. P. J. Chem. Phys. 2006, 124, 
044103. 

Constantin, L. A.; Fabiano, E.; Della Sala, F. Phys. Rev. 
B 2011, 84, 233103. 

Fabiano, E.; Constantin, L. A.; Della Sala, F. J. Chem. 
Phys. 2011, 134, 194112. 

Buhl, M.; Kabrede, H. J. Chem. Theory Comput. 2006, 
2, 1282-1290. 

Zhao, Y.; Truhlar, D. G. J. Chem. Theory Comput. 2005, 
1, 415-432. 

Fabiano, E.; Constantin, L. A.; Della Sala, F. J. Chem. 
Theory Comput. 2014, 10, 3151-3162. 

Jurecka, P.; Sponer, J.; Cerny, J.; Hobza, P. Phys. Chem. 
Chem. Phys. 2006, 8, 1985-1993. 

Marshall, M. S.; Burns, L. A.; Sherrill, C. D. J. Chem. 
Phys. 2011, 135, 194102. 

Luo, S.; Zhao, Y.; Truhlar, D. G. Phys. Chem. Chem. 
Phys. 2011, 13, 13683-13689. 

Peverati, R.; Truhlar, D. G. J. Chem. Theory Comput. 
2012, 8, 2310-2319. 

TURBOMOLE, V6.3; TURBOMOLE GmbH: 
Karlsruhe, Germany, 2011. Available from 
http://www.turbomole.com (accessed August 2014). 
Weigend, FFurche, F.: Ahlrichs, R. J. Chem. Phys. 
2003, 119, 12753. 

Weigend, F.; Ahlrichs, R. Phys. Chem. Chem. Phys. 

2005, 7, 3297-3305. 

Paier, J.; Marsman, MKresse, G. J. Chem. Phys. 2007, 
127, 024103. 

Constantin L. A.; Pitarke, J. M. Phys. Rev. B 2011, 83, 
075116. 

Horowitz, C. M.; Constantin, L. A.; Proetto, C. R.; 
Pitarke, J. M. Phys. Rev. B 2009, 80, 235101. 

Grimme, S.; Antony, J.; Ehrlich, S.;. Krieg, H. J. Chem. 
Phys. 2010, 132, 154104. 

Tao, J.; Perdew, J. P.; Staroverov, V. N.; Scuseria, G. E. 
Phys. Rev. Lett 2003, 91, 146401. 

Perdew. J. P.; Ruzsinsky, A.; Csonka, G. I.; Constantin, 
L. A.; Sun, J. Phys. Rev. Lett. 2009, 103, 026403-026407; 
(E) Phys. Rev. Lett. 2011, 106, 179902. 

Ruzsinszky, A.; Sun, J.; Xiao, B.; Csonka, G. I. J. Chem. 
Theory. Comput. 2012, 8, 2078-2087. 
del Campo, J. M. ; Gazquez, J. L. ; Trickey, S. B. ; Vela, 
A. Chem. Phys. Lett. 2012, 543, 179. 


92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 

111 

112 

113 

114 

115 



